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We use the AdS/CFT conjecture to investigate the thermalization of large- N c M = 4 Super Yang- 
Mills plasma in the limit of large but finite 't Hooft coupling. On the gravity side, we supplement the 
type IIB supergravity action by the full set of C(a' 3 ) operators, which enables us to derive 0(A _3//2 ) 
corrections to the emission spectrum of prompt photons in one model of holographic thermalization. 
Decreasing the coupling strength from the A = 00 limit, we observe a qualitative change in the way 
the photon spectral density approaches its thermal limit as a function of the photon energy. We 
interpret this behavior as a sign of the thermalization pattern of the plasma shifting from top/down 
towards bottom/up. 



Introduction. Figuring out the precise way, in which 
quark gluon plasma created in a heavy ion collision ap- 
proaches local thermal equilibrium prior to its eventual 
freeze-out and hadronization, is a question of utmost im- 
portance for a successful description of the system. In 
the limit of weak coupling, the thermalization pattern 
is known to be of the bottom/up type, with the soft 
excitations reaching the thermal limit before the hard 
ones PQ; recently, it has been found that the driving 
force behind this behavior is related to plasma instabil- 
ities [21 E]. In contrast, studies of thermalization in the 
strongly coupled limit of conformal field theories, utiliz- 
ing the gauge/gravity duality [1HS], have pointed towards 
top/down type thermalization [7Hl2j (for an overview of 
holographic calculations in the heavy ion context, see also 
|13j ) . This indicates the probable existence of a transi- 
tion between the two behaviors at intermediate coupling. 

The range of 't Hooft couplings relevant for a heavy ion 
collision, A = g 2 N c ~ 20, falls clearly outside the realm 
of weak coupling techniques (see e.g. j!4l for a discussion 
of this issue). Starting from the opposite end, A = 00, 
one might however wonder, if a consistent inclusion of 
strong coupling corrections might allow one to access this 
region. In the case of large- N c J\f = 4 Super- Yang-Mills 
(SYM) theory, these corrections are available through the 
inclusion of the full set of 0(a' 3 ) type IIB string theory 
operators in the supergravity action (151 116) . resulting 
in corrections to different physical quantities typically 
proportional to the parameter 7 = |C(3)A~ 3 / 2 [T7H20"] . 
These calculations are often conceptually straightforward 
but technically very involved, one prominent example be- 
ing the recent determination of the (thermal) photon pro- 
duction rate beyond its well-known A = 00 limit [2T] in 
refs. [251 [23]. In these articles, the authors first took ad- 
vantage of their earlier results [24] to derive 0(A~ 3 / 2 ) 
corrections to the Schrodinger equation of a bulk U(l) 
vector field, and then proceeded to solve this equation 
to obtain the photon spectral density on the field theory 
side. 

In an altogether different generalization of ref . [21] , two 
of the present authors studied photon and dilepton pro- 
duction in an out-of-equilibrium J\f = 4 SYM plasma 



[251 l2"tj] , working within a model of holographic thermal- 
ization that involves the gravitational collapse of a thin 
spherical shell in AdSs space [271128]. This model is ho- 
mogeneous in three-space and allows an analytic treat- 
ment of photon and dilepton spectral densities in the 
limits of high frequencies and late times. In these two 
cases, characterized by the slow relative motion of the 
shell, one may namely approach the dynamics of the sys- 
tem in a quasistatic approximation, where the shell is 
treated as a static object when formulating the bound- 
ary conditions of different bulk fields. This presents an 
important simplification in the calculations, and in par- 
ticular allows one to verify the validity of the fluctuation 
dissipation theorem, necessary to connect the field theory 
spectral densities to the corresponding (photon/dilepton) 
production rates. 

In the paper at hand, our plan is to combine the ap- 
proaches of refs. [22J [231 US] to study the production of 
prompt photons in a thermalizing M = 4 SYM plasma 
beyond the usual A = 00 limit. In particular, we aim to 
extract the relative deviation of the photon spectral den- 
sity from its thermal limit in a specific out-of-equilibrium 
state, and to examine the dependence of this quantity on 
photon energy at different couplings. As we will discuss 
in the following sections, we argue that this information 
can be used to study the pattern, with which plasma 
constituents of different energy approach local thermal 
equilibrium. 

The setup. We work within a model of holographic 
thermalization that involves the gravitational collapse of 
a thin shell of matter in AdSs space. The radial co- 
ordinate of the shell (in a coordinate system where the 
boundary is located at r = 00) is denoted by r s . It satis- 
fies at all (finite) times the relation r s > r^, where r/j is 
the Schwarzschild radius of the eventual black hole, cor- 
responding to the temperature of the field theory in its 
final thermal state. The dynamics of the shell, i.e. the 
dependence of r s on the field theory time t, can in princi- 
ple be solved from the corresponding equation of motion; 
cf. the discussion in [25] . For the purposes of the present 
paper, we however leave this question aside, as we will 
in any case work within the quasistatic approximation, 



2 



where the thermalization process is most conveniently 
parametrized by the values of r s /r/j. 

Inside the shell, the metric of the spacetime is that of 
pure AdSs , while outside it is given by an AdS black hole 
solution; the effects of being out of thermal equilibrium 
then become manifested through the changing boundary 
conditions of different bulk fields at the location of the 
shell. The model was initially introduced in ref. [27] and 
further analyzed e.g. in [TJ][55] (see also [55] for a discus- 
sion of how the collapsing shell can be realized by turning 
on a scalar source on the boundary). Our notation and 
conventions follow those of ref. [25] , to which we refer the 
interested reader for more details. 

In order to study photon production, we next add 
to the JV = 4 SYM theory a gauge field coupled to a 
conserved current corresponding to one of the U(l) sub- 
groups of the SU(4) R-symmetry of the theory. As dis- 
cussed in [21] . this is most conveniently performed in a 
way, in which all Af = 4 scalars and fermions obtain equal 
charges, leading to an anomaly free theory. To leading 
order in the electromagnetic coupling a and to all orders 
in A, the photon production rate per unit volume is then 
given by the formula 
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where 115^ is the electromagnetic current Wightman 
function. In thermal equilibrium as well as in the qua- 
sistatic limit of the falling shell setup (see the discussion 
in appendix B of [IS]), the fluctuation dissipation the- 
orem further allows one to relate H< to the transverse 
photon spectral density Xfiv> 

»rn< (*° = k) = n B (k°M(ko) = n B (k°)x(ko) , (2) 
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This implies that it suf- 



fices to consider the retarded correlator of the electro- 
magnetic (in practice, R symmetry) current in the SYM 
theory, which simplifies the dual gravity calculation in a 
dramatic way. 

As we want to study photon production to the next- 
to-leading order in a small curvature expansion, we must 
include C(a' 3 ) corrections to the type IIB supergravity 
action. This results in [l5l [16] 
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where Rio is the 10-dimensional Ricci scalar, cf> the dila- 
ton, and F§ the five-form field strength, and where we 
have set the curvature radius of AdSs space to unity. In 
the 0(7) terms, C denotes the 10-dimensional Weyl ten- 
sor, while T stands for a complicated tensor built from 
the five- form and defined e.g. in eq. (5) of |22) . 

The 7-corrected AdS black hole metric derived from 



the above action has the form [T51 [T71 [TB] 
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where f(u) = 1 — u 2 and u = r^/r 2 is a dimensionless 
coordinate, in which the boundary of the AdS space is 
located at u = 0. The different functions appearing here 
are given by 

K{U) = e 7 [»(«)+4&(«)l j = &1 K«) , L(u) = e~< ^ , 
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while the 7-corrected relation between and the field 
theory temperature reads r/, = nT/(l + ^yf^)- 

The calculation. To determine the required R current 
correlator including all necessary 7-corrections is in prin- 
ciple a very laborious exercise. Fortunately, we can use 
here the results of refs. [22] [24J [30] , in which the authors 
consider a vector field living in the above background 
metric. Specializing to the momentum space component 
of a transverse electric field, Et(u,w), with u the cor- 
responding frequency, they arrive at a Schrodinger-type 
action 
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where the boundary term reads $(u) = ^'(u)^(u). The 
relation between "J and Et can be found from [30] and 
has the form *(tt) ee Z(u)E t (u), S(u)- 1 ee l/y/f(u) + 
jp(u), with the function p(u) given by 
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With these definitions, the equation of motion for the 
field takes the simple form *"(u) - V(u)^(u) = 0, 
where the potential V(u) reads 



V(u) = 
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Working in the collapsing shell model, we need to solve 
the above equation for * outside the shell, i.e. for u < 
u s . Inside the shell, u > u s , neither the (pure AdS) 
metric nor the equation of motion for ^ is altered by 7- 
corrections [31] [35] , and thus the solution to the latter 
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FIG. 1: The photoemission rate dF-y/ko, normalized by 
aNcT 3 and evaluated for r s /rh = 1.01 and A = 
oo, 120, 80, 40 (from bottom to top at the peak value). The 
qualitative form of the curves stays unaltered for all small and 
moderate values of r s /rv 



can be read off e.g. from ref. [55]. Here, one only needs 
to recall the relation between the frequencies measured 
inside and outside the shell, 



f m = f(u s )K 2 (u s ), (9) 



where the — subscript refers to the inside and + to the 
outside. The two solutions furthermore need to satisfy a 
set of junction conditions at the shell [25] . 



= VU^+(u s )/^{u s ) , 
*'_(«.) = / TO w (tt + (u)/E(tO)| 
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which result from demanding continuity of the physical 
electric field and where the it-derivative on the left hand 
side of the second relation does not operate on the pa- 
rameter u s in w_ . 

Following the above steps, we first obtain for the inside 
solution 



*_(u) = (u) + (u) , 



(11) 



where can be read off from eq. (10) of [55J (noting 
that inside the shell ^(u) — Et(u)) and where the cor- 
rection term is available through a simple expansion 
of the modified Bessel functions. For the outside solution, 
we on the other hand write 



-(*? n + 7*L) + c+(*L+7^ 
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where the subscripts 'in' and 'out' refer to infalling and 
outgoing field modes at u = 1. The zeroth order solutions 
can again be read off from ref. [55] (after a rescaling by 
V f( u ))i while the correction terms are obtained through 
a numerical solution of the equation of motion discussed 
in the previous section, subject to the relevant boundary 
conditions. Matching this outside solution to the inside 
one via eq. ( 10 1 then provides us with a 7-corrected result 



for the ratio c_/c + , through which we have determined 
the behavior of ^(u) in the entire AdS space. 
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FIG. 2: The photoemission rate dTj/ko, normalized by 
aN^T 3 and shown for A = 100 and r s /r h = 1.1, 1.01, 1.001, 1 
(from bottom to top at very small frequencies). 



Results. To obtain the photon spectral density \ 011 
the field theory side, the above results simply need to be 
inserted into the relation I3"3l 
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which we expand to linear order in 7. From here, we 
easily obtain both the relative deviation of the spectral 
density from its equilibrium limit, 



R(u) 



X(&) ~ Xth(w) 
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as well as the corresponding photon production rate, 
dF-v ak 
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Our results for these quantities are to be compared on one 
hand to the A = 00 off-equilibrium results of ref. [26] and 
on the other hand to the 0(a' 3 ) corrected equilibrium 
calculation of ref. [55] . 

In fig. [l] we first display the photoemission rate for 
different values of A (chosen such that the strong cou- 
pling expansion is still applicable), with the shell always 
residing at r s /rh = 1.01. The pattern we observe is very 
similar to that described in the thermal case in ref. [22j : 
Decreasing the 't Hooft coupling from A = 00, the peak of 
the spectrum increases and moves towards smaller to. At 
the same time, fixing the value of the coupling to A = 100 
and varying r s /rh is seen to largely reproduce the quali- 
tative findings reported in ref. [53] for A = 00; see fig. [2] 
above. The main difference between the behavior of our 
curves and those plotted in fig. 3 of [55] is the faster shift- 
ing of the peak of the finite-A spectral density towards 
larger w, when Tg/r^ is increased. This has the effect of 
suppressing the value of the peak due to the appearance 
of the Bose-Einstein distribution in eq. ([T]), best seen in 
the blue r s /rh = 1.1 curve in fig. [2} 

In fig. [3] we display our perhaps most prominent find- 
ing, the behavior of the relative deviation of the spectral 
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FIG. 3: The relative deviation R, plotted for r s /r h = 1.01 and A = oo, 500, 350 (left) and A = 150, 100, 75 (right), with the 
amplitudes of the curves increasing with decreasing coupling. While in the A = oo case the amplitude of the oscillations gets 
damped at large u>, for all finite values of A it first decreases but ultimately starts growing linearly with u>. 



density from its thermal limit, R(6j), for r s /r^ = 1.01 
and A = oo, 500, and 350 (left), as well as 150, 100 and 
75 (right). The quantity is observed to exhibit oscilla- 
tions similar to those discussed already in [26], which 
increase in frequency and decrease in amplitude as the 
shell approaches the horizon. Interestingly, the behavior 
of the fluctuations as functions of oj seems to depend on 
the value of A rather strongly: As soon as one leaves the 
strict A = oo limit, the asymptotic large-w behavior of 
the amplitude of R changes from a 1/u suppression to 
a linear enhancement. For a given value of A, the am- 
plitude of R is in fact observed to have a minimum at 
some oj — w m i n (A), signifying the frequency, for which 
the spectral density is closest to its thermal limit. This 
is in direct contrast with the usual geometric picture of 
strong coupling top/down thermalization, where one ex- 
pects the hardest modes to always thermalize first. A 
numerical study further shows that within the validity 
of the strong coupling expansion, the function u; m i n (A) 
always obeys a power law behavior w m i n (A)/T ps cvA* 3 , 
where the positive constants a and /3 depend on the value 
of r s /rh only very mildly. For r s /rh — 1.01, we obtain 
a ps 1.25, /3 ps 0.69. 

Although it is good to recall that the photon frequency 
oj does not directly correspond to the energy of the 
plasma constituents, it is tempting to speculate that the 
above behavior of R might be indicative of the thermal- 
ization pattern of the plasma changing from top/down to- 
wards bottom/up as the value of A is decreased. After all, 
it is at the very least clear that the highly energetic pho- 
tons with oj ^> T cannot have been emitted by the soft 
excitations of the plasma. Furthermore, one should note 
that our observation regarding the shift in the asymp- 
totic behavior of R is extremely robust. First, it is seen 
at all non-infinite values of A, including a regime where 
the strong coupling expansion is guaranteed to converge. 
And second, going beyond the limit of a static shell would 
result in corrections proportional to 1/(toj), where r is 
the characteristic time scale related to the motion of the 
shell. As discussed in section 4 of [25] , these contributions 



are expected to be strongly suppressed for large enough 
values of oj, such as those considered in fig. [3] above. 

Conclusions. In the paper at hand, we have used the 
AdS/CFT conjecture to study the behavior of an out-of- 
equilibrium large- N c AT = 4 SYM plasma at large, yet 
finite 't Hooft coupling. The calculation was carried out 
within one particular model of holographic thermaliza- 
tion, where the approach of the plasma towards thermal 
equilibrium is modeled via the gravitational collapse of 
a thin shell of matter in AdSs space. In this setup, we 
were able to derive the leading 0(1/A 3 / 2 ) corrections to 
the photon production rate and the corresponding spec- 
tral density. Our most interesting finding, displayed in 
fig. [3] was interpreted to reflect a change in the thermal- 
ization pattern of the plasma from top/down towards 
bottom/up as A is decreased. In particular, our results 
indicate that for values of A of relevance to real life heavy 
ion collisions, one may already be quite far from a strict 
top/down behavior. 

Although our results offer only indirect evidence of the 
pattern, with which strongly coupled SYM plasma ther- 
malizes, we feel that they constitute an interesting start- 
ing point for further study. An immediate generalization 
of the present calculation should be a similar determina- 
tion of energy momentum tensor correlators within the 
collapsing shell model, providing direct information on 
the behavior of the plasma constituents themselves. On 
top of this, it would be intriguing to extend the present 
calculation beyond the quasistatic approximation, as well 
as to other, more realistic models of holographic thermal- 
ization. These studies would be crucial to assess the uni- 
versality of our findings, as well as the conclusions to be 
drawn. If it turned out that the observed change in the 
thermalization pattern of a non-Abelian large- N c plasma 
with decreasing A was a generic prediction of holographic 
calculations, it would clearly highlight the importance of 
including strong coupling corrections to the present stud- 
ies. At the same time, such a result would be highly 
encouraging for the future prospects of modeling ther- 
malization in heavy ion collisions using holography. 
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